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TECHNICAL NOTE 5018 

A THEORETICAL STUDY OF THE EFFECT OF FORWARD SPEED ON THE 
FREE-SPACE SOUND -PRESSURE FIELD AROUND PROPELLERS 
By I. E. Garrick, and C. E. Watkins 

SUMMARY 


The -so-und-pressiire field of a rotating propeller in forward fll^t 
in free space is analyzed hy replacing the normal-pressure distribution 
over the propeller associated with thrust and torque by a distribution 
of acoustic pressure doublets acting at the propeller disk and subject 
to uniform rectilinear motion. The basic element used to synthesize the 
field is the pressure field of a concentrated force moving uniformly at 
subsonic speeds, for which an expression generalizing one of Horace Lamb's 
for the fixed concentrated force is given. This result is presented both 
for the moving and fqr the fixed observer. The strength of the doublet 
distribution is related to the thi*ust and torque distribution in a con- 
venient blit approximate way. The sound field is expressed by integration 
over the propeller disk, and also by integration over an effective ring, 
and is given both for the near pressure field and, in a simpler form, 
for the far field. Known results for the zero-forward-speed case present 
themselves in the special case of Mach number M = 0. Some illustrative 
exaii^iles are calculated and discussed. 


INTRODUCTION 


The rotating propeller is the source of an Intense sound-pressure 
field which can be associated with the periodic reactions on the medium 
■arising from the distribution of pressure rotating along with, the blades. 
This pressure distribution consists in part of a distribution due to 
thickness of the blades, whose resultant force in subsonic potential flow 
is zero, and in part a distribution due to angle of attack and camber 
of the blades whose integrated effect Includes the Induced drag and corre- 
sponds almost wholly to the thrust and torque distribution over the blade. 
Another source of propeller noise may be associated with flow separation 
and with friction or shear due to the boundary layer; both effects lead 
to vortlcity shed into the wake and hence the designation vortex noise. 

The vortex noise and the noise due to thickness (where wave drag is not 
a large factor) are, however, for actual propellers normally of a consid- 
erably smaller magnitude than the rotational sound due to torque and 
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thrust} hence only the latter effect will he considered in the present 
work. 


A large number of investigators have studied various phases of the 
determination of the sound or noise field of rotating propellers. In 
addition to the references cited, a bihliography is Included of repre- 
sentative work on this subject. A singplification that has frequently 
been made in propeller-noise investigations is to limit the consider- 
ations to the static or standing propeller. Ihe work of Gutin in 195 ^ 

(ref. 1) represents a development^of this type that makes possible a 
satisfactory prediction of the amplitude of sound pressure due to thrust 
and torque of a propeller rotating on a stand in still air. Althou^ 

Gutin 's theory is applicable to the near oscillatlng-pressure field of 
the propeller, his results in reference 1 are limited to the determi- 
nation of the fimdamental and the first few harmonics at a distance far 
from the propeller, that is, several diameters away. The determination 
of the near pressure field, however, has been of concern both frcm struc- 
tuiral and physiological considerations. Hubbard and Regier (ref. 2 ) 
extended the application of Gutin *s work to describe the oscillating 
pressure field and to determine the ait5>litude of noise at points near 
the propeller, in some cases within a blale chord length from the tip. 

They Investigated analytically the effect of several of the parameters 
that enter in the theory and also gave con^jarisons with experiment which 
were quite satisfactory. 

The existing theoretical work has found useful application for "I 

static conditions and for conditions of low forward speed, for example, 
near take-off. For conditions of high forward speeds, however, many 
pertinent questions have arisen as to the possible effects of the forward 
speed on the oscillating pressure field, of the propeller. A few investi- 
gators have examined phases of this problem} Bryan, Hart, Shirokov, 

Blokhintzev, and particularly Kuessner and Billings (see bibliography) 
may be mentioned, but most of this work seems incomplete or difficult to 
apply. It appealed desirable, especially for applications, to reexamine 
the theoretical problem for the propeller at forward speed so as to have 
it arise as a straightforward generalization of existing work for the 
zero -forward- speed case. The purpose of the present paper is therefore 
to extend the theory of Gutin (ref. 1 ) and the work of Hubbard and 
Regier (ref. 2 ) so as to include the effect of subsonic forward speed 
of the propeller on the near and the far oscillating sound-pressure 
field, caused by torque and thrust of the propeller. 

This paper Includes the following material: (a) Hie sound- pressure 

field associated with a uniformly moving concentrated force is given with 
details of the explicit" development contained in the appendix. This 
result serves as a basic element for synthesizing the total pressure 
field. It is expressed both for the case of the observer or field point 
considered to be moving along rectilinear ly and uniformly with the pro- • 

peller in free space and for the case of the observer considered fixed 
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and the propeller in uniform flight, (h) llie disturhance forces associ- 
ated with the pressure distribution acting on the medium in the plane of 
the propeller disk are presented in the manner similar to that of Gutin. 
(c) The so\md -pressure field resulting from the combined use of (a) and 
(h ) is shown in the form requiring integrations over the propeller disk 
and in the simpler form making \iae of an effective propeller radius. 
Approximations valid for the far field, which reduce to Gutin 's results 
for the case of zero forward speed, are also given, (d) Some numerical 
examples are calculated and discussed. (e ) Remarks on the usefulness 
and limitations of the analysis. 


SYMBOIS 

A(r) chordwlse distribution of thrust acting on a radial element 

of a propeller blade 

B number of propeller blades 

b(r) width of propeller blades 

c velocity of sound 

Cp power coefficient, p/pn^D^ 

Cqi thrust coefficient, T/pn^D^ 

D propeller diameter 

P(r) chordwlse distribution of forces perpendicular to the thrust 

of a propeller blade and giving rise to the torque 

Fg conponents of force vector F 

Bessel function of first kind with index mB 
Mach number, V/ c 
order of harmonic 
propeller rotational speed, rps 


k = “ 




M 

m 

n 



k 
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p 

power 


f* 

p 

pressure 



lpl 

pressirre magnitude 




pressure due to torque 



Pj 

pressiire due to thrust 



p 

■^rms 

root-mean -square pressure 



Q 

torque 



E 

length of propeller blades 



Re 

effective length of propeller "blades 



r 

radius to a "blade element 



r, 9 

polar coordinates in yz-plane 


A 

s = 

(x - xj_f + ^2[(y - y^f + (z - zi)^ 



So = \ 

fx2 + 32 y 2 



^ = ]f( 

- ^if + (y - yi)^ + (z - zif 



So = \ 

(x2 + y2 



So^ ® 

polar coordinates in xy-plane 



T 

thrust 



t 

time 



V 

forward velocity 



X, y. 

z Cartesian coordinates 






fm. 
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M. 


p = ^/l - 


Pr 

blade angle 

T1 

propeller efficiency 

P 

fluid density 
\ . 


M(x - + S 


psrlod, 2rt/Bft 

angular velocity 

frequency of mth harmonic, mBfi 

fundamental frequency, Bft 


T = 


rfl 


a 

(a 

“1 


ANALYSIS 

Extension of a Formula in Lamb's "Hydrodynamics" for the 
Sound-Pressure Field of a Fixed Concentrated Force 
to That of a Moving Force 


On the basis of acoustic considerations of the classical hydrodynamic 
equations. Lamb (ref. 5) gives the pressure at any field point x, y, z 

associated with an external periodic force Ae^“^ acting in the x-direc- 
tion and concentrated at the location x^, y^^, z-|_ as 


icu(t-f) 


P = T— 


A _3 

8x s 


where 


= \|(x - xif + (y - yi)^ + (z - zi)' 


( 1 ) 


( 2 ) 
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The concentrated force may he regarded as stemming from the limit of a 
distribution of an increasing pressure difference over a decreasing 
area, whose product in the limit of zero area is equal to the force. 

For an arbitrary concentrated force of components Fy, F^ the result 

generalizes as Lamb indicates to 

8 - 1) ^ a - 1) ^ a - f) | , , 

Sbc s 5y 8 8z 8 



This formula shows that the pressure field obtained from the concentrated 

•1 

force has the character of a doublet or a dipole j in its periodic form 
it has played a central role in the development of the noise and pressure 
fields of propellers due to the propeller torque and thrust. However, 
the formula refers to a concentrated force and to coordinates fixed in 
space. To treat the case involving effects of forward speed, it is expe- 
dient to obtain the required extension of the preceding formula. 

The extension of eqioation ( 5 ) required for the case of a uniformly 
moving concentrated force is given in equation (4), which follows. 

Details for its derivation are supplied in the appendix. Let the con- * 

centrated force have components F,,., Fy, Fgj let it be moving uniformly 

with velocity V in the direction of the positive x-axis; and let the a 

coordinate system also move imlformly with the same velocity. Then, 
with X, y, z now denoting the field point referred to this coordinate 
system and with x^_, y^^, Zj_ denoting the coordinates of the concentrated 

moving force — 


P 


hit 


:(^ - f) 


8 x 


a - f) , 8 - I) 

c»y S 8z S 


(h) 


where 


S = ^(x - xi)2 + p2|^y _ y^^2 + 
M^x - Xq_) + S 

P = VI - 


(5) 




Since the acoustical pressure is p = -p the dipole is an accel- 

ot 

eratlon doublet or the time rate of change of a flow doublet. 
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Figure 1 Illustrates the geometric meaning of the quantities s, S, 
and a. lEhe force located at and the field point at 

Q(x,y, z) are hoth i n un iform motion vlth velocity V In the x-dlrec- 
tlon. The distance OQ = s. The Influence at Q at time t stems 
from the action of the force vhen It was at position P, where P Is 
obtained from the relation OP/OQ = M, the forward-speed Mach number. 

Ihe distance PQ gives then directly the phase radl\is cr. The perpen- 
dicular dropped from 0 onto PQ determines QR, which Is equal to the 
amplitude radius S. For M = 0 both S and a reduce to the ordinary 
radius s. Con^jarlson of equations ( 5 ) and (i(-) shows that the effect of 
the forward speed leads to replacement of s by S In the amplitude and 
by a In the phase. 

The force variation of main concern herein Is that of a harmonically 
periodic force F having con 5 )onents F^^., Fy, F^ varying In time as e^^. 
Equation (4 ) becomes 


P = 






( 6 ) 


where k = m/c. Equation ( 6 ) is the basic relation to be made use of in 
the subsequent analysis. 


Disturbance Forces In the Propeller Plane 

Consider as does Gutin (ref. l) a propeller disk oriented so that 
the axis passing throu^ the center of the disk coincides with the x-axis 
and let the propeller be ass\jmed to rotate, as In figure 2 , In the 
yz-plaue (x = O) with positive values of x corresponding to points 
ahead of (and negative values of x corresponding to points behind) the 
propeller disk. The propeller Is considered to move uniformly with 
velocity V In the positive x-direction. Points In the propeller disk 
are designated by j polar coordinates (as shown in 

fig. 2 ), by 


y^ = r cos 0 
z^ = r sin 0 






( 7 ) 


For definiteness, let the propeller be rotating counterclockwise as seen 
by an observer looking into the propeller toward the slipstream. 
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Each element of the propeller is acted on by the surface pressure 
distribution and this distribution may be resolved into a thorust force 
in the direction of the axis of rotation (the x-axis ) and into a force 
associated with the torque which acts about the axis of rotation opposing 
the rotation. Equal and opposite reaction forces to these are exerted on 
the medium. The points of application of these forces are imagined to 
act in a single plane designated as the plane of rotation or as the pro- 
peller disk. 

The reaction of the surface presstire distribution of the rotating 
propeller on the medi^Im at instant is to be replaced by fixed 
periodic forces acting at the propeller disk. The propeller disk itself 
will be considered covered with the necessary singularities in the pres- 
svire or acoustic radiators of proper strength and harmonic content to 
correspond to this normal-pressure distribution. These singularities 
will be seen to be acceleration sources (acoustic radiator of zero order ) 
for those forces which act symmetrically on both sides of the blade, that 
is, whose net force over each element of the blade is zero, and to be 
acceleration doublets (acoustic radiator of order one ) for those forces 
acting antisymmetrically, that is, whose net force over each element of 
the blade corresponds to the difference in pressure over both sides of 
the blade. The doublet distribution is that which is needed to repre- 
sent the thrust and torque distribution, in particular that part of the 
thrust and torque distribution associated with pressures acting normal 
to the blade surface. This part is practically all of it, except that 
arising from skin friction. (The effect of blade thickness may be taken 
into account by introducing either flow sources or acceleration sources: 
the flow sources would lead to the sound-pressure field due to thickness 
by the classical "piston" effect of the moving blade on the flow field; 
the acceleration sources would deal directly with the assumed or known 
contribution to the pressure distribution over the blade due to thickness.) 

Consider an element of the propeller at distance r from the axis; 
let dr be its radial length and b its width measured in the projec- 
tion onto the plane of rotation. Let the forces acting on the propeller 
element on each blade be A(r)dr in the axial flight direction and 
F(r)dr in the direction opposed to the direction of rotation. Equal 
reaction forces acting opposite to the fli^t direction and in the direc- 
tion of rotation, respectively, are exerted on the medium. These quanti- 
ties are related to thrust T and torque Q by the relations 


^ = BA(r) 
dr 


^ = BrF(r ) 
dr 


( 8 a) 
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pR 

T = B / A(r)dr 
'-I 0 

pE 

Q = B / rF(r)dr 


where B Is the number of "blades and R the length of each "blade. 

The periodic impulse or reaction experienced at any element of the 
disk may be expanded in a Fourier series. For simplicity let the element 
considered first be located on the radial line 0 = 0 (the y-axis ) and 
afterwards be located arbitrarily. To start with, assimne that the forces 
are uniformly distributed over the projected width b of the blade 
element, that is, that the distribution of pressure difference over the 
blade chord is rectangular. Then in the area element r dr d0 of the 


plane of rotation, the forces A(r )dr and F(r )dr 

b b 


act on the 


medl'um during the time interval in which this element is eclipsed by the 
projection of the propeller element. If the overlapping of the element 

starts at t = 0, it will end at t = t = — and the overlapping of the 

element will start again by the next blade at t = Tq = — , where n is 

the angular velocity of the propeller. The rectangular -type forces 
experienced at the element of the disk located at 0=0 by its periodic 
eclipse by the blade may be developed in a Fourier series: 


Fi*(t) 


A(r ) ^ dr for 0 < t < t 

b 


for T < t < Tr 


and similarly 




Fp*(t) 


Fp*(t) 


cos(mBfit - em) 

m=l 


F(r) dr for 0 < t < t 

b 


for T < t < Tf 


Bq + ^ ®m cos(mBat - 
n^l 


( 10 ) 
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where 


= ~ A(r) — sin mtt — dr d0 

“*■ nut - 


B_ = — F(r) — sin nut — dr dB ^ 
^ DM h ( 




T _ mBflT 
To " 2 Tq 


( 11 ) 


The constants Aq and Bq which correspond to the instantaneous average 

thrust and torque over the "blade element and to the associated momentum 
shed into the slipstream do not give rise to sound and need not he 
expressed. The phase angles and which are small for the lower 

harmonics, are needed to preserve generality of discussion since Aj^ cos 
corresponds to the coefficients of pure cosine terms and Aj^ sin ejjj to 

coefficients of pure sine tenns in the Povirier series. For the assumed 
rectangiilar distribution, and % may he observed to "be equal to 

zero if the origin of time is chosen at the overlapping of the center 
line of the blade by the y-axis since then only a pure cosine series 
suffices to express the distribution. 

The general formulation of interest, where the disk element r dr d0 
is located at an angle 9, may be expressed directly from these results, 
for, on a second area element r dr d0 shifted with respect to the first 
by the angle 0 in the rotational direction, there act periodic forces of 
the same magnitude, but retarded by the time 0/n. The corresponding 
Fomrier developments are 


( V 00 “ 

t - ^ } = Aq + cos(mBflt - mB0 - em) 

m=l 

FgCt) = F 2 *(t _ = Bq + cos(mBnt - mB6 - T] ) 


( 12 ) 


The quantity 


(15) 
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is small for the lower harmonics, especially for hlade elements not near 
the hub. Blade elements near the hub, where r is small, are eliminated 
from consideration since they contribute very little to the air forces. 

One may then replace in equation (U) sin by go that, noting 

To Tq 

the relations in equations ( 8 a), 


Ajji — A(r )dr d9 

It 

B, « ® F(r)dr d9 

^ 3t 


dT 

« — dr d 9 
dr 


» i ^ dr d 9 
r dr 




Uk) 


This approximation, which wei^ts the harmonic content of all the har- 
monics equally, may be observed to correspond to the case of the thrust 
and torque distributed over a zero blade width (that is, to the mathe- 
matical pulse sometimes termed the Dirac delta function). This approxi- 
mation becomes relatively less valid when exceeds about ^ 

To ^ 

^sin ^ = 0 . 707 } ^ = 0.785^ or when hm. ~ exceeds about unity; that Is, 
when the order of the harmonic exceeds — or, roughly ^ l/solldlty. 


Moreover, the assimiptlon of a uniform rectangular distribution of 
the forces across the blade has been made for convenience, and other 
distributions may be treated if desired. Some of the possible errors in 
the assumption leading to equation (l^l-) ane discussed briefly by Gutin 
in reference 1 and are shown to be generally negligible for the lower 
harmonics. Regier and Hubbard (ref. i^) also discuss this assun^jtion in 
an illuminating manner. Figure 3 ^ which is taken from reference h, 
shows the relative hannonic content of different assumed distributions 
of the same total load for; (a) the mathematical sharp pulse, (b) a 

triangular "hat" pulse = O.03, (c) a rectangular pulse -i- = 0 . 03, 

2 rtr 2 Ttr 

and, (d) a rectangular pulse = 0 . 06. Blade widths of actual pro- 

2.TCT 

pellers tend to be between cases (c) and (d) for the most effective parts 
of the propeller disk. The sharp-pulse assumption generally tends to 
overestimate the magnitude of the higher harmonics. When the proper 
distribution is known, appropriate correction factors for the required 
harmonic may, of coirrse, be applied to the magnitude of the results given 
by the pulse solution. 
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The formulation in equation ( 12 ) need not he limited to the one in 
which the propeller force distribution is uniform thro\igd‘°^'t the pro- 
peller cycle. If interest should be attached to a nonuniform distribu- 
tion, as occurs for the propeller yawed or pitched with respect to the 
flight path or for one experiencing interference effects, it is readily 
possible to allow for these effects by permitting the amplitude of the 
distribution to become a function of 9. 


The Sound-Presstire Field 

The general expressions for the resolution of the forces associated 
with the thrust and torque on a radial blade element (for example, 
eqs. ( 12 )) may be put in the usual convenient complex form (whose real 
part may correspond to the formulation of interest): 



(15) 


where cd = mBft and where for the special case of the rectangular pulse 
distribution the coefficients are given as in equation (ll). 

With the use of the approximations for Aj^ and given in equa- 

tion (l4), the periodic forces acting on the mediimi at the element r dr d9 
of the disk can be expressed for any given harmonic in the forms 


jt dr 


rt r dr 


dr de 


( 16 ) 


T, 1 1 dQ , i(aJt-mB0), 

F„ „ —2; cos 0 e dr d0 

jt r dr 


In the derivation of eqiiatlon (I 6 ), the phase angles ej^ and T]j^, which 

can generally be neglected for the lower harmonics, have been put eqioal 
to zero, a value which corresponds in the case considered (as has been 
mentioned) to the choice of 0=0 to correspond to the overlapping of 
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the center line of one of the blades at the time ’ t = 0 . Ihe index m 
indicates the order of the harmonic considered. As discussed in the 
preceding section, appropriate factors depending on the harmonic nim- 
her m, functions of the radius, and, for unsymmetrical loading, fiinc- 
tlons also of 0 may (if known) he applied to the terms in equation ( 16 ). 


Pressure relations involving integrations over the propeller disk . - 
The pressure at any point in the free-space field produced hy these com- 
ponents of the uniformly moving periodic forces F^, Fy, F^ is given in 

equation ( 6 ). With the use of the components of force given in equa- 
tion ( 16 ) and hy integration over the propeller disk, the total oscil- 
lating pressiire p for any given harmonic m (the index m will he 
hereinafter dropped) is obtained from equation ( 6 ) as a sum of the pres- 
sure due to thrust p^ and that due to torque Pq as follows: 

P = + Pq ( 17 ) 


where 


1 r® ^ i(mt-mB 0 ) 


A. 

2bc S 


dr d 0 


( 18 ) 




^ ^ gi(cot-mB0) 1 

0 J 0 rdr® rSe 


,-ika 


dr d 0 


( 19 ) 


and 

for 

The 

out 


where, in the expression for p^^, the operation 
convenience to replace the equivalent operation 


1 A 

r 30 
sin 0 


has heen used 


— COS 0 

oy oz 


indicated differentiations with respect to x and 0 can he carried 
hy use of the expressions for a and S given in equations ( 5 ) 


3 / IkM ikx x \ 

3x S S p 2 q2 J 


1 3 e-i^^ 


e-lka 


S 



(y sin 0 - 


( 20 ) 


r 30 S 


S 


z cos 0 ) 
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Hence, 


P 


T “ 



dl _imB0 e~^^^/lkM 
dr S [^2 


Ikx 

• 


— \dr de 


and 


( 21 ) 



1 ^ 
r dr 


-lmB9 

S U 



sin 0 - z cos 6 )dr cL9 


( 22 ) 

Another somewhat siit 5 )ler expression for p_, can be obtained from eq.ua- 
tion ( 19 ) by Integration by parts with respect to 0 


P 


Q 


e^^ r R r 1 ^ -imB0 e~^^^ 
'J 0 0 ^ ^ ® 


imB dr d0 


(25) 


Pressure relations involving effective ring approximation . - Appreci- 
able simpllficatlpn may be achieved for calculation purposes by making 
use of the approximations inherent in the assumption of an effective pro- 
peller radius Rg so that the integration with respect to r is avoided. 

Equations (21) and (2j) then reduce to 


and 


glcut y^2rt 


^/ikM 

\p2 0 % gpj S 





g-1 (mB0+ka) 
S 


d0 


(24) 


( 25 ) 


where, in both S and a, the points y^ and z^ have the values 

= Re cos 0, Zj_ = Rg sin 0. In effect, equations (24) sind ( 25 ) imply 
that the propeller disk has been replaced by an annular ring in which 
the entire thrust and torque are concentrated. The effective radius of 
this ring varies somewhat, as calculation shows, with the load distri- 
bution and with the order of the harmonic. Deming (ref. 5) has shown by 
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calculation in special cases for the static propeller that the ring 
approximation is a reasonably good one. An effective radius of the order 
of 0. 8 r is considered reasonable for adaption in Initial calculations. 

The magnitude of the root -mean-square pressure ^ is of interest, 

since most sound-recording Instruments are calibrated in terms of it. The 
contribution of the harmonic of order m to the root-mean-square pres- 
sure is 


P 


rms 




or 


■ rms 


8it^ 


(26) 


where 


p 2jc 

r 

— 

— 1 


L ' 

— — cos(mB9 + ka) + 

L" 

7 " I) 

_ qEL. 

sin(niB0 + ka)» 


d9 

S 


B = 




— 

— 1 


L 1 

sin(mB9 + ka) + 

T hfM + 2'\ 
_ 32V sj 

Bell 

cos(mB0 + ka)- 


S 


When rotational symmetry exists as for the condition of symmetrical 
loading, it is convenient in making the numerical calciilatlons to let 
the field point be in the xy-plane, an arrangement which can always be 
attained by suitable choice of the line 0=0. 

Pressure relations- for the far field . - A fvirther simplification in 
the results can be reached if the distance s from the propeller disk 
(0,yi, to the field point (x,y,0) is large. Then from 


s 


2 


+ (y - yj_f + 



there is obtained 


o 


yyi 


s 


o 


s s 


(27) 
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where Sq = yx^+y^ (the dlsteince from the field point to the center 
of the propeller dish). It also follows from 

+-p^(y - yi)^ + 

that 

pStt 

S « Sq . -—i: (28) 

C>o 

where 

So = \/x2 + p2y2 


If terms of order i/s^ are neglected in compaxison with terras of 
order l/S, the sound pressure due to thrust is obtained from eqvia- 
tion ( 2 ^)-) by use of equations (5) and (28); 


iojt -i^(So+Mx) p 2 ii _imB0+ifc|-EeCOs0 


From the known Integral relation 

^2TC 


nelK 

/ de = 2rti%(X) 


it follows that 


Pt 


= T l“® e 

2rtS„ 


Similarly, from equation (25), 


( 350 ) 


9(“ • tW") 


Pq = 


,mBfl _ ^ /kyR, 


Q ^^mB.1 p, 

Re^ 


JmE 


( 352 ) 
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Observe that the argument of the Bessel function may he replaced hy 


kyRe 


- "®Mrot I- 


(53) 


where = —^} the L&ch number corresponding to the rotational speed 

at the effective radius. (The far -field approximation and the effective- 
radius approximation made use of in arriving at eq.s. ( 51 ) and (52) need 
not, of course, he made simultaneously, since the required integration 
of the Bessel functions with respect to the radius may he carried out, 
at least numerically, without difficulty. ) 

Introduce the fundamental frequency denoted hy = BJi, so that 


o) = mBJ2 = kc = mo)]_ 


The pressure magnitude for any harmonic m is then given hy 




|P 



(3^) 


This result may he compared with that for the case of zero forward speed 
(M = O) given in reference 1; 



mo)]_ 

2rtcsQ 



- Q 


Be 

(OiRe^ 



(35) 


It can he noted that the forward-speed Mach number 
torque term containing Q only in the replacing of Sq 
distance from the field point to the propeller huh) hy the smaller dis- 
tance Sq = \/x2 + p2y2^ This substitution occurs both in the argument 
of the Bessel function and in the outside factor 1 /Sq. The thrust term 

containing T is more strongly affected, since, in addition to this 
change, the field effect ahead of and behind the propeller disk is Influ- 
enced hy the term M + — corresponding to a backward shift hy = -M 

bo 

Increase because of the factor 



af fects th e 
= + y^ (the 


and also hy an 
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APPLICATION TO A SPECIFIC PROPELLER 


In order to give some indication of the effect of forward-speed 
Mach number on the sound pressure of a propeller, calculations based on 
equation ( 26 ) for the near field and on equation ( 5 U.) for the far field 
have been made . For this purpose, a two-blade research propeller having 
a 10 -foot diameter and operating under the various conditions summarized 
in table I was chosen. As may be noted from this table, the propeller 
is assumed to operate at constant power; that is, the-power coeffi- 
cient Cp and the torque Q are held constant. As the forward-speed 
Mach number is changed, the blade angle and consequently the 

thrust T are changed according to the propeller charts of reference 6 
so as to make the assumed conditions consistent with actual test oper- 
ating conditions. Only the fundamental, the first harmonic m = 1, is 
considered in these examples and the value chosen for the effective 
radius Rg is 0.8R or 4- feet. It should be pointed out that the sound 
pressure computed from the data of table I would be obtained in pounds 
per square foot, since these data are given in English imlts. These 
values have been converted to dynes per square centimeter in the results 
given in the figures by multiplication by the approximate conversion 
factor 480. 


Calculations for the near field . - Calculations of the root-mean- 
square pressure = -^|p|^ based, on equation ( 26 ) are made for 


various values of x in the range from -O. 5 D to 0 . 5 D along the line 
y = 0.6 d (6 feet), that is, along a line l/5 of the propeller radius from 
the tip and extending a distance of 1 radius behind to 1 radius ahead of 
the plane of rotation. Hiese results are showi plotted as a function of 
x/D for several Mach numbers in figure 4. It can be noted in this fig- 
ure that, for each Mach number, two peaks generaJly appear, one ahead of 
and one behind the plane of rotation. For the set— of conditions under 
consideration (see table l), the hipest peak amplitudes of pressure 
occur just behind the propeller plane at" values of x/D in the range 
from -O.I 5 to -0.075* As the Mach number is increased from 0 to 0.4, 
the peak amplitudes decrease in magnitiide, but as the Mach number is 
increased from 0.4 to 0.9 this trend is reversed and the peak amplitudes 
increase in magnitude, the peak for M = O .9 being about 1.4 times that 
for M = 0. Also, as M increases from 0 to 0 . 9 , the point at which 
the highest peak pressure occurs moves somewhat nearer the propeller 
plane. Thus, the generally severe soxand -pressure conditions at take-off 
(M R* 0 ) tend to be alleviated in flight at the lower Mach numbers but may 
be reached again and even exceeded at the higher flight Mach numbers. 

The calculations for low Mach numbers presented in figure 4, it should 
be noted, are in substantial agreement with both calculated and measured 
results of reference 2. The hi^ peak pressures obtained for the highest 
Mach numbers indicate that, with propeller -driven airplanes, the sound 
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pressures generated near the tips of the propellers are of significance 
and of possible concern with regard to both structiiral considerations 
and passenger comfort. 

Calculations for the far field. - Calculations of the root-mean- 
sq.'uare pressirres based on equation (5^) for various values of x along 
the line y = 2D (20 feet) are shown plotted in figure 5- trend with 

regard to Mach nimber in this figure is about the same as was noted in 
figure 4 for y = 0. 6 d, but the relative effect of the forward Mach num- 
ber appears greater at a distance. However, because of the greater dis- 
tances in figure 5 than in figure 4, the peak an 5 )lltudes are considerably 
less. 


Directional characteristics of the sound field . - For some purposes, 
especially with regard to calciilations for the far field, it is desirable 
to consider the character of the sound pressmre in terms of polar coordi- 
nates. For this purpose, substitutions can be made in eqiiation (5^) 

(and eq. ( 2 . 6 )) as follows: 


X = Sq cos <6 y = Sq sin $ s^ = yx^ + y^ 


( 56 ) 


QIhe quantities x/Sq and v/Sq in eqtiatlon ( 3 ^) (and eq. ( 26 )) 
may then be written as 


X 

COB $ 

So 

^1 - M^sln^J 

y _ 

sin $ 

So 

\/l - ^^sin^$ 


With these substitutions, equation (3^) becomes 


(37) 


P = 


2 rtcs 


moj]_ 

- m 2 ; 


s I 


f» . ^ 

^1 - M^sin^t^ 


- Q 


Be 




kRg sin O 


^mB 




M'^sln'^®/ 


( 38 ) 

Calculations based on equation ( 38 ) for the Mach numbers of 0 and 
0.8 for a constant value of Sq = 2D are shown plotted as dashed curves 
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in figure 6 . The solid cuinres in this figure represent results of calcu- 
lations along the line y = 2D obtained by replacing s^ by — - — 

^ sin $ 

with y = 2D in equation (50 ) for the various Mach numbers considered 
in figures 4 and 5 - A con^iarison of results of calculations along the 
line y = 2D with those along the circle Sq = 2D for M = 0 and 

M = 0 . 8 indicates that"; for_htg^ Mach numbers (of the order of M = 0. 8 ), 
the peak pressures calculated along a line y = Constant-' are about the 
same as those calculated along the circle with radius equal to the con- 
stant value of y. Observe the second pressure peak which has developed 
in the forward location at M = 0. 9« 

Sep arate ccaigionents due to torque and to thrust of the sound field . - 
In order to give some indication as to the nature and proportion of 
sound pressure associated with each of the quantities, thrust T and 
torque Q, the root-mean-square pressures ' associated with_ each of these 
quantities for M = 0 and M = 0.8 are plotted for y = 0.6 d (6 feet) in 
figure 7* These plots show, as do equations ( 26 ) and (5^), that the root- 
mean-square presstires associated with the propeller torque are symmetri- 
cally distributed with respect to the plane of the propeller for all Mach 
numbers, whereas the root-mean-square pressures associated with the 
thrust are symmetrical with respect to this plane only for M = 0. For 
the particular propeller and operating conditions imder consideration, 
the amplitudes of pressure associated with thrust for low Mach numbers 
are higher than those associated with torque, but- for hign Mach nvimbers, 
the opposite is true. In the Interpretation of figure J, it should be 
recalled that the results depend on the assumed operating conditions and 
that the torque Q is the same ( 2,680 ft-lb) in both parts of the figure; 
whereas the thrust is 1,850 pounds at M = 0 and 510 pounds at M = 0.8. 
Plots of this type can be used to obtain the sound pressure for various 
thrust and torque coefficients for a given propeller, since these coef- 
ficients appear as factors in equations ( 26 ) and ( 5 ^) aJ^d hence can be 
normalized. 

It is of particular interest to note in the plot for M = 0.8 that 
the peak pressiires associated with both thinst and torque are- consider- 
ably greater than the peak pressure associated with the sum of these 
quantities. This result indicates that' the phase relationship between 
the two congjonents is, in this case, such that each has a canceling 
effect on the other. Also in this case, the pressure associated with 
thrust may have its greatest" value ahead of ■the plane of the propeller; 
whereas the total pressure has its greatest" value behind the plane of 
the propeller. 

The discussion of. the figures based on the specific examples illu- 
strates the fact that the results do depend markedly on the assumed 
operating conditions. Moreover, only the fundamental harmonic (m = l) 
for the specific two -blade propeller has been Illustrated. Trend studies 
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on effects of higher harmonics, number of hlad.es, and d.lfferent oper- 
ating cond-itions would he of considerable interest. Some preliminary 
calculations on the soimd pressures associated with the fundamental of 
a six-hlade propeller have shown a greater relative effect of the 
forward-speed Mach number. 


CONCLUDING REMARKS 


Expressions have heen given for the sound -pressiure field due to the 
distribution of thrust and torque for any given harmonic of a rotating 
propeller in miiform subsonic flight. The general expressions (eqs. (IT ) 
to (25 ) ) involve integrations over the propeller disk; approximate 
expressions for the near field (eqs. (2^-) to (26)) and for the far field 
(eqs. (51 (52), and ( 3 ^)) involve integrations over a ring with an 
effective radius. The numerical exan5>les have illustrated some free- 
space sound -pres sure results for the fundamental of a specific two-blade 
propeller \mder various operating conditions at various forward speeds. 

It is pertinent to remark again on some of the limitations of the 
analysis. The analysis presented utilizes the torque and the thrust 
distributions, which may be given empirically or theoretically, in such 
a way as to require that they arise purely from pressures acting normal 
to the surface of the blades. Empirical values of the thrust and torque 
Include a contribution, generally small, due to skin friction and to 
separated flow. Some caution is then needed in the use of the results. 
For example, the empirical tor-que term is Bomewhat larger than the torque 
due p\irely to normal press'ures and, hence, the associated sound result 
due to torque may be overestimated sli^tly; similarly, the empirical 
thruat tezm may slightly iinderestlmate the sound due to thrust. At high 
tip speeds there may be significant contributions to the noise due to 
wave drag associated with the thickness. These contributions are mainly 
taken into account by the effects of the wave drag on the torque. (Other 
soxmd effects of the thickness should be separately calculated and 
Included but these effects have not been explicitly presented herein. ) 

Another way of looking at the approximation, and perhaps a generally 
more convenient one for the study of trends, is to consider that the 
assumed torque and thrust distributions are actual theoretical ones, 
obtained if necessaiy by adjustment of the blades of the propeller in a 
potential flow; and hence the sornid -pressure field is that corresponding 
to the chosen thriist and torque, but in potential flow. Although the 
soimd -pressTzre field and the aerodynamic velocity field have been con- 
sidered as separate, it is of Interest and of significance that the same 
concepts leading to the calculation of the soimd field (made use of in 
the form of the acceleration potential) can lead to the linearized aero- 
dynamics of the propeller in compressible flow, including the represen- 
tation of the vortlcity left behind in the wake as a result of spanwise 
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variation of loading. The theoretical Induced drag and the theoretical 
wave drag (which is a form of acoustic loss) are i^erently included in 
the representation employed. 

It may he worthy of repetition that the pressure formulas as given 
tend to overestimate the contribution of the hi^er harmonics and that 
appropriate factors based on chordwise loading can be devised and may 
be required. Moreover, because of the theoretical (as well as the 
enplrical) change in aerodynamic loading along the blade radially with 
increasing tip Mach niunber and forward-speed Mach number, the appropriate 
factors will change. In addition, the effective radius will be altered 
with the loading. Thus, althou^ 0.8 blade radius may be suitable as an 
effective radius for normal loading, a smaller value for the effective 
radius may be more suitable for conditions where unloading of the tip 
occ\irs. 

A calculation study of trends under different assumed conditions, 
the effect of unsymmetrical loading, or of dual propellers, effects of 
obstacles or boundaries on the free-space resiilts, and experimental 
confirmations for the In-fli^t propeller are interesting matters for 
further investigation. 


Langley Aeronautical laboratory. 

National Advisory Committee for 
Langley Field, Va. , August 


Aeronautics , 

51, 1955. 
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APPENDIX 

THE SOUND-PRESSURE PIEED ASSOCIATED WITH 
A UNIFORMLY MOVING FORCE 
Derivation of Equations (il-) and (6) 


Lamb (ref. 3 ^ P- 502 -) derives from hydrodynamical principles tbe 
differential equation satisfied by the press\ire condensation s* in the 
acoustic field associated with arbitrarily vaiying fixed forces acting 
on the medium 


8y^ 



= + M.+ 

Sy hz) 


(Al) 


where Tamb 'a X, Y, Z are actually the extraneous forces per unit volume 
divided by the density p. (in effect. Lamb is dealing here with the 

acceleration potential. ) In terms of the perturbation press\are p = pc^S'* 
the equation is 

C^p = div F 




5y 8z 


(A 2 ) 


where F is the arbitrary force per unit volume having coniponents 

F^. Lamb shows that if the periodic force Ae^^ is imagined 
concentrated on an infinitely small space at (x]_,y2_, and to be in the 

direction of x, the pressiare at the field point (x,y,z) (the distance s 
from the location of the force) is given by 


P 



^ s 


(A 3 ) 


where 


s . ^(x - xi)2 + (y - y^)2 + (z _ z^)2 
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so that the concentrated force is equivalent to a double source or 
acceleration doublet whose- axis is in the direction of the force. For 
a general harmonically periodic force F having conponents I‘’x^ ^yt 

varying in time as e^^^ 


P = 




s 




8 \ 

bzj s 


(A4) 


where k = — . 

c 


For an arbitrary time -dependent concentrated force F(t) located 
"the pressure at field point (x,y,z) can be expressed as 


P = 5- dlv 


- I) 


1_ 

4n 


A ~ i) + 1. ~ t) ^ ~ I) 

8x S ^ s 8z B 


(A5) 


where the differentiations affect the result only through the variable s. 

The extension to a moving concentrated force may, within the frame- 
work of small-perturbation theory, be made in several ways - for example, 
as mentioned in reference 7* One formal procedure utilized by Kiissner, 
for example, makes use of the invariance properties of the wave equation 
and utilizes combined Galilean and Lorentz transformations. The following 
direct procedure given briefly in reference 7 is believed to be of 
intrinsic interest;^ as it represents a simplification in Prandtl's pro- 
cedirre (ref. 8) for the case of moving constant- source distributions, a 
procedure which consists of schediiling a succession of fixed sources in 
a path to act consecutively one sifter the other so as to represent in 
effect the desired so-urce moving along its path. Let the arbitrary con- 
centrated force act only as an impulse diiring an inf inj-beslmal in-terval 
at time t = T. The impulse may be -written as 


(a6) 


F(t)6(t - t) 
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where the inrpiilse f\inction S(t) = 0 for t ^ 0 and is characterized 
for T = 0 as having unit area with respect to t. The useful property 
of the Impulse function of "sifting or selecting" a value of a function 
Is exhibited hy the following relation (see, for example, ref. 9 } P* 61): 



/ CO 

f (r) 5 (a 


T)dT 


= f(cT) 


(AT) 


Ijet a succession of such Impulses act, one following the other, in a path, 
points along which are given In space-fixed coordinates by 


t = I(t) ti = ti(t) 5 = ^(t) 

The effect at time t of all such impulses which act before the time t 
is then given from. Lamb ' s result (eq., (A 5 ) ) as 



where, as has been defined, there is a nonzero contribution to the 
integral only for values of t defined by the characteristic relation 


t - T = |- = ^||{x - |(t)]^ + [y - + [z - (A9) 

which expresses the distance between the source point and the field point 
in terms of the time of travel of the outgoing waves. 

The integral in eqmtlon (A8) corresponds to a summation of temporary 
fixed sources. To represent the case of uniform rectilinear motion with 
velocity V in the positive x-dlrection, let the so-urces be located on 
the l-axis and flow consecutively one after the other at the positions 


| = Vt tj =0 S = 0 (AlO) 

so that for t = -« the source is located at | = -00 and for t = 0 
the source is located at the origin. It follows that the distance 
between source point and field point is 


s = ^(x - Vt)^ + y^ + 


(All) 
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Replace in equation (A8) the variable t by 0 vhere 



(A12) 


By the sifting property of the 5-function (eq. (A7 ) ) the integral will 
have a value only for 0=0 and equation (A8) becomes 


P = 


7^- dlv 


F(t) ^ 
s d0 


0=0 


(AI 5 ) 


where the q\iantlty within the bracket is now to be evaluated. Equa- 
tions (All) and (A12) give a quadratic equation for r which, on choice 
of the solution that leads to t < t and with 0=0, gives 



where 

S = ^(x - Vt)^ + p^(y^ + z^) 

p = )ll - 



(Al4) 


Equations (All) and (A12) also give 


1 ^ _ c 

s d0 cs - V(x - Vt) 


where replacing s by its characteristic value c(t - t) from equa- 
tion (A 9 ) and replacing t by its value given in equation (Ali)-) yields, 
for 0 = 0 , 


s d0 S 


(AI 5 ) 


where S is defined in equation (Aik). Equation (AI 5 ) becomes 


F 



V_ S\_l^ 
c2 "" " cj^2 


(AI 6 ) 


p = i dlv 


s 
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This result Is referred to space-fixed coordinatesj it is finally- 
desired also to express this result in terms of a field point (xq^Yq, Zq) 

of a coordinate system moving uniformly along with the source located at 
yj_, Zj_ so that one may put 


X - Vt = Xq - X]_ y = Yq - y^ Z = Zq - ZjL 


After this substitution is made, the zero subscripts may be dropped to 
yield as the end result corresponding to equation (A 5 ) for the source 
and the coordinate system in imiform rectilinear motion; 


P = 



(ait) 


where 


S = \j{x - Xi )2 + p 2 |^y . ^ (2 - z^)^ 

M^x - x^) + S 


P = \/l - 


Eq-uatlon (AI7 ) corresponds to the result given in equation (i<- ) of 

the analysis. For the periodic harmonic force F varying as e^^, it 
may be expressed in the form 

p = i (F - V) (AI8) 

which is equivalent to equation ( 6 ) of the analysis. Equations (AI7) 
and (A 18 ) are the sou^t-for generalizations of equations (A5) and (a 4 ) 
of the appendix, for the case where the disturbance and the field point 
are in -uniform rectilinear motion. Con^jarlson of these equations shows 
at once that s is replaced by S in the amplitude and by a in the 
phase. A geometric interpretation of these quantities is shown in fig- 
ure 1 and discussed in the analysis following eqixatlon (5). 
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Remarks On the Case of the Moving Disturbance 
and the Fixed Observer 

A few remarks are in order on the significance of equation (A 16 ). 

In this equation, the field point fx,y, z) is given in space-fixed coordi- 
nates while the disturbance force (or propeller) is in the in-flight 
condition moving with velocity V in the positive x-directlon. Its 
location is given by Xj_ = Vt, y^ = 0 , z^_ = 0 so that at — t =0 it is 

located at the origin. It may be observed that ilie distance between 
disturbance and field point S is numerically the same whether given 
by the relation in eqmtion (AI6) or (AI7). Hence, the pressure magni- 
tude at any observer location for a so\ind -radiating element of the 
in-fll^t propeller is llie same as that for the observer moving along 
with the propeller, provided the proper instantaneous distance between 
observer and propeller is used. However, there will be a difference in 
the frequency perceived by the observer. This frequency will be that 
for the case of the uniformly moving observer modified by the Doppler 
effect; thus, the frequency of- each harmonic is modified by c/(c + 
where Vj- is the conqjonent of the propeller -forward speed (with proper 

sign; minus, for approaching, plus for receding from, the observer) 

In the acoustical direction from the observer to, say, for the far 
approximation, the hub of the propeller. The acoustical direction is 
not quite that from the observer to the propeller location; it actually 
points from the observer to the location of the propeller when the sound 
which reaches the observer was emitted (direction GEP rather than QO 
in fig. 1). The Doppler frequency factor 

m ^ 

OP 


c^^c -I- Vj.) is given geomet- 
S 
a 


rlcally in figure 1 by the ratio of 
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TABLE I 

SEA- LEVEL OPERATING CONDITIONS FOR A lO-FOOT-DIAMETER 
TNO-BLADE PROPELLER (DATA FROM REFERENCE 6 ) 

[p = 815 hpj Cp = 0.10) Q - 2,680 Ib-ftj k = 0.2968^ 


M 

V/nD 

deg 


■ n , 

percent 

ft -lb 

0 

0 

24 

0.11 

0 

1,850 

.1 

.42 

26 

.11 

46 

1,850 

.2 

.84 

28 

.094 

75 

1,600 

.5 

1.26 

53 

.065 

85 

1,090 

.k 

1.68 

37 

.C 52 

87 

875 

.5 

2.10 

42 

.041 

86 

690 

.6 

2.52’ 

hi 

.052 

80 

580 

.7 

2.94 

50 

.025 

70 

420 

.8 

5.36 

55 

.0185 

62 

510 

.9 

3.78 

56 

.0122 

46 

205 
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Flg\ire 1.- Geometric representation of e, S, and a 






Relative amplitude of harmonics 
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(a) Sharp pulse. 


(h) Triangular pulse. = 0.03 • 

23tr 


(c) Rectangular pulse. — ^ = 0.03.- 

23tr 


(d) Rectangular pulse. = 0.06. 

2 jtr 

Figure 3 .- Effect of inpulse shape on the relative amplitudes 

of the harmonics. 
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Flg\Jire 4.- Distribution of the root -mean-square pressures for the 
fimdamental (m = l) of a tvo-hlade 10 -foot- diameter propeller at 
several forward-speed Mach numbers. Cp = 0.10; k = 0.29686; 
y = 6 feet. (Operating conditions of propeller given in table I.) 
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Figure 5-- Distribution of the root-mean-sq.\iare pressures for the 
fundamental (m = l) of a two-hlade 10-foot- diameter propeller at 
several forward-speed Mach numbers. Cp = O.lOj k = 0.29686} 
y = 20 feet. (Operating conditions of propeller given in table I.) 
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(a) M = Oj T = 1850 Ib. 


(b) M 0.8; T = 5IO lb. 


Figure 7‘- Distribution of the root-mean-sq.uare total presBuree and the root-mean-aquare pres- 
sures associated vith thrust and torque for a two-blade, 10-foot-diameter propeller at forward 
speed Mach numbers of 0 and 0.8. Cp = 0.10; m = 1; k = 0.29686; y = 6 feet, (operating con- 
ditions of propeller given in table I.) 
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